PHASE SPACE MODELS FOR STOCHASTIC NONLINEAR PARABOLIC 
WAVES: WAVE SPREAD AND SINGULARITY 



ALBERT C. FANNJIANG 

Abstract. We derive several kinetic equations to model the large scale, low Fresnel number be- 
havior of the nonlinear Schrodinger (NLS) equation with a rapidly fluctuating random potential. 
There are three types of kinetic equations the longitudinal, the transverse and the longitudinal with 
friction. For these nonlinear kinetic equations we address two problems: the rate of dispersion and 
the singularity formation. 

For the problem of dispersion, we show that the kinetic equations of the longitudinal type 
produce the cubic-in-time law, that the transverse type produce the quadratic-in-time law and that 
the one with friction produces the linear-in-time law for the variance prior to any singularity. 

For the problem of singularity, we show that the singularity and blow-up conditions in the trans- 
verse case remain the same as those for the homogeneous NLS equation with critical or supercritical 
self-focusing nonlinearity, but they have changed in the longitudinal case and in the frictional case 
due to the evolution of the Hamiltonian. 



1. Introduction 

In this paper we consider nonlinear Schrodinger (NLS) equation with a random potential 
(1) z— ^(z,x) + ^Ax^'(z,x)+7-ig|^'|2-^(z,x)+^F(2L„xL,)^(z,x) = 0, 

X G M"', a > 

where 7 is the Fresnel number (defined below), is the linear coupling coefficient for the random 
potential V, which is rescaled by two large parameters and Lx, and g is nonlinear coupling 
coefficient with g > representing the self-focusing case and g < the self-defocusing case. Here a is 
a positive constant and a = 1 corresponds to the cubic NLS equation. We are particularly interested 
in the regime of low Fresnel number 7 <C 1 with a rapidly oscillating potential Lx,Lz, n ^ 1. 

Finite-time blow-up or wave collapse is a well-known effect for the self-focusing, (super)-critical 
{da > 2) NLS equation without a random potential jSH]. In this case the nonlinear focusing 
effect dominates over the linear diffraction effect of lA. The question remains whether the added 
random potential and its scattering effect would prevent the formation of singularity. Given the 
wide range of scales present in this problem the numerical simulation as well as theoretical analysis 
is undoubtedly extremely challenging. Although current numerical results, e.g. jS], have already 
indicated that a white-noise-in-z potential without the feature of rapid fluctuation in x does not 
prevent blow-up. Indeed, without the self-averaging effect of a x-rapidly fluctuating potential the 
large z-fluctuations in the white-noise potential may drive the system to a state with low, negative 
Hamiltonian, thus developing singularities in very short time. These singularities in our view are 
the small scale, randomly fluctuating singularities which may be delayed or removed in the presence 
of a x-rapidly fluctuating potential. 

To capture the robust, large scale blow-ups we derive several phase-space model equations cor- 
responding to different scalings (in 7, /i, L^;, L^). We use these kinetic equations to elucidate the 
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problems of singularity formation and rate of dispersion. The main ingredient of our analysis is 
the variance identity for these phase space transport equations. In general, we show that random 
scattering changes the blow-up conditions for the homogeneous NLS eq. but does not prevent 
singularity formation on the large scales (see Conclusion for the discussion and summary of our 
results). We also give various upper bounds for the blow-up time. 

As a simple by-product of the variance identity, we also derive various bounds on the rate of 
dispersion for random Schrodinger waves, including exact expressions in the critical case, da = 2. 
We show that depending on the random scattering mechanism, one can have the cubic (z^), the 
quadratic (z^) and the linear (z) rate of dispersion prior to singularity, if there is any. To our 
knowledge, these are significant improvements over current results (e.g. jH], [21], [Z]) which are 
mostly for the linear problem. 

The cubic NLS equation is a model equation in nonlinear optics, describing the spatial distribu- 
tion of the stationary electromagnetic field in a nonlinear medium with the Kerr effect |2Uj . |15j . 
Here let us give a brief derivation of the NLS equation with a random potential. 

The electric field E with a fixed polarization in a lossless medium satisfies the Helmholtz equation 

V|S(x) + k^r?E{i) = 0, X E 

where k the wavenumber and n(x, \E\'^) > is the refractive index of the medium. For a weakly 
nonlinear and weakly fluctuating inhomogeneous medium we write 

n(x, \E\'^) = n(l + n(x)) + n2\E\^ 

with n <C 1 and the Kerr coefficient n2 being small. Consider the approximate equation 

(2) V|£;(x) + k^n^Ei^l) = -2n^n{SC)E{Sl) - 2nn2\E\^E{Sl). 
The electric field is taken as 

£'(z,x) = ^'(z,x)exp (inkz), x = (2;,x) G 

where z is the coordinate in the direction of propagation and x is the transverse coordinates. The 
modulation ^ is assumed to be slowly varying on the scale k~^. 

Substituting the ansatz into eq. Q and making the parabolic approximation we obtain for the 
modulation ^ the equation 

d 1 

(3) ikn—'if{z, x) + - ^/^^(z, x) + k'^nn2\^\'^^{z, x) + k'^n^n{z, yi)^{z, x) = 0. 
Let us write 

n{z, x) = vV{z, x) 

where <C 1 is the standard deviation of n and V{z, x) is a zero-mean, z-stationary, x-homogeneous 
random field. Let Lx,Lz 3> 1 be the two length scales of the wave beam in the transverse and 
longitudinal directions, respectively and 

x = x/L^., z = z/L;,, 

be the rescaled variables. Let us introduce the dimensionless parameters 

H = knLzV, 7 = 7^^, g = L:,n2k-f. 
kriL^ 

Writing the wave field ^ in the new variables x, z, dropping the tilde and rescaling we obtain 
with a = 1. 

The critical case of the cubic NLS with d = 2 and g > Q is our primary example; another is the 
sub-critical case of a planar waveguide (thin film) with d = 1. The opposite case of a self-defocusing 
Kerr nonlinearity occurs in semiconductor waveguide jl9j . |18j . Similar equations also arise in many 
other contexts such as the Langmuir wave in plasma and Bose-Einstein condensation, see |23] and 
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references therein. Therefore, we wih formulate the results for the general power-law NLS eq. 
with transverse dimension d > 1. Because the variable z is time- like we will refer to it as "time" 
occasionally, especially in the discussion of finite-time singularity. 

2. WiGNER DISTRIBUTION AND WiGNER-MOYAL EQUATION 

We consider several families of scaling limits, as 7 ^ 0, L^, L^^ ^ oo, which are first distinguished 
by whether 

(4) e = -fL^^o or e = 1. 

This, of course, is not sufficient to ensure the existence of scaling limit until we specify the strength 

of fj.. 

Our phase-space model equations for the low Fresnel number regime is based on the Wigner 
equations. We consider the Wigner distribution of the form 

(5) ^(x, p) = ^ I e-P-y* [z, X + ^) ^* (z, X - ^) dy. 

The Wigner distribution has many important properties. For instance, it is real and its p-integral 
is the modulus square of the wave function, 

(6) / Ty^(x,p)dp = |M'(x)|2^p(x), 

SO we may think of l^(x, p) as wave number-resolved mass density. Consequently, the mean 
and variance Vx oi x are given by, respectively. 



xVK(x, p)dxdp, 
K = j |x-xpVF(x,p)dx(ip 



where Sx is the second moment of x 



Sx = j |xpVF(x,p)dxdp. 



Also, we have that 



/ W^(x,p)dx = (^)'^|$-|2(p/^). 

and the energy flux 

(7) 1(^V^* -^'*V^') = / pW%x,p)dp 

2« jRd 

SO that the mean and variance of p, Vp, are given by, respectively 

p = J pW{x,p)dxdp, 



where Sp is the second moment of p 



Sp = j \p\'^W{x,p)dxdp. 



3 



In view of these properties it is tempting to think of the Wigner distribution as a phase-space 
probabihty density, which is unfortunately not the case, since it is not everywhere non-negative (it 
is always real- valued though). 

It is straightforward to derive the closed-form equation for the Wigner distribution [H] 

dW 

(8) — + p-V^W+U-yW + nVW = 0, 

oz 

with the Moyal operators 

(9) U^W{x,p) = iy"e*'i-"7"MW^(x,P + 7q/2)-VF(x,p-7q/2)][/(z,q)dq, U = g\^\^'' 

(10) VPF(z,x,p) = if^J e*''"''^ [VF,(x,p + 0q/2) - H^,(x,p - 0q/2)] ^(zL^dq) 

where U, V denote the Fourier transforms of U, V, respectively. Formally we see that as 7 ^ 

(11) U^W{x, p) ^ UoW{x, p) = V^U{z, x) • Vpt^(x, p). 
We use the following definition of the Fourier transform and inversion: 

•^-'^(x) = I e'^-^g{p)dp. 

We shall refer to eq. (jSJ as the nonlinear Wigner-Moyal (NWM) equation. We only need to consider 
the weak formulation 

(12) ^ J eWdj^dp - j V- VxOI^dxdp - j U^QWdyidp - V@Wdy.dp = 

for smooth, rapidly decaying test functions 0. 

One advantage of working with eq. is that one can use it to evolve the mixed-state initial 
condition, instead of the pure-state one given in Q. This is important in the context of modeling 
quantum open systems. A mixed state Wigner distribution is a convex combination of the pure- 
state Wigner distributions © described as follows. 

Let {^'q} be a family of functions parametrized by a which is weighted by a probability 
measure P{da) . Denote the pure-state Wigner distribution Q by VF[^']. A mixed-state Wigner 
distribution is given by 



(13) J W[^a]P{da) 



The limits as 7 —> of the mixed state Wigner distributions constitute the so called Wigner 
measures which are always positive |12| . jl3j . |17j . Evolution by eq. ()12() preserves the form 
In particular, for such initial data we have 

(14) jw{z,x,p)dp > 0, VxGM^ 

(15) jw{z,x,p)dx > 0, VpGM'^ 

Multiplying (|12|) by W and integrating by parts we also see that the evolution preserves the 
L2(M2rf).norm of W, i.e. 

— / \W\'^dxdp = 0. 
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Let us be more explicit about the random potential V{z,x.). We assume that z-stationary, 
x-homogeneous random field V is square-integrable and admits the spectral representation 

V{z,x) = J exp {ip ■ x)V{z,dp) 

with the z-stationary spectral measure V{z, •) satisfying 

E[V{z, dp)V{z, dq)] = d{p + q)$o(p)dpdq. 

Here E denotes the ensemble expectation and the function ^oIp) is the transverse power spectrum 
density which is assumed to be rapidly decaying as |p| or |p| — > oo. The transverse power 
spectrum density is related to the full power spectrum density ^(w,p) in the following way 

^o(p) = j '^{w,p)dw 

with <I'(k) satisfying <I>(k) = <!>(— k),Vk = (tf,p) £ R"'^"'^. Without loss of generality we will also 
assume 

(16) ^{w, p) = ^{-w, p) = ^{w, -p) = ^{-w, -p). 



3. Phase space transport equations 

3.1. Linear kinetic equations. First we summarize what has been established in the linear case 
when U {z, x) is a given function. In this case we refer to eq. © with a given U as the linear 
Wigner-Moyal (LWM) equation. The principal feature of the scaling is the separation of scales in 
the given potential U and the random potential V as they appear in the LWM equation, c.f. @ 
and (|lUj) . Under pretty general conditions (the most important of which being the integrability of 
the maximum correlation coefficient of V{z, •) as x-homogeneous-field-valued z-stationary process) 
one can prove that, as 7 ^ and /i, L^, 00, the weak solution of the LWN equation converges 

in law to the weak solution of the linear Boltzmann (LB) equation or the linear Fokker-Planck 
(LFP) equation, described below, depending on whether 9 also tends to zero or not [H],^n] (see 
also 0). 

Linear Boltzmann equation {6 = 1): 



d 

(17) a^^^-^' ^' P) + P • ^xM^(^, X, p) + Vxt/ {z, x) • VpW{z, X, p) 

= 27r y" K{p, q) [Wiz, x, q) - W{z, x, p)]dq 

with a nonnegative kernel K(p, q) given respectively as follows. 

(a) : If ~ \/T^, Lx <C Lz then 

(18) K(p,q) = cl>(0,q-p) 

(b) : If /i ~ /L^, Lz-^Lx<. d > 3 then 



(19) i^(p,q)=<5( 
(c): If /i ~ ^/Ll, Lx ~ Lz then 



|qp-|pp. 



<I>(w, q — p)dw 



(20) 



K(p,q) = cI>(l3LJPl!,q_p). 



(23) 



Linear Fokker-Planck equation (9 —>^ 0): 

(21) ^W{z, X, p) + p • V^W{z, X, p) + V^U{z, x) • VpW{z, x, p) 
oz 

= Vp-DVpTF(z,x,p) 

with a symmetric, nonnegative-definite matrix D given as follows. 

(a) : If /i ~ 6~^y/L^, <C then 

(22) D = Try $(0,q)q(g)qdq. 

(b) : If ^ ~ 9-^y/lZ, Lz < < > 3 then 

D(p) = 7r|p|"W ^{w,p±)dw p±0p_i_dp_i_. 

J p px=0 J 

(c) : If /i ~ 9^^ y/L^, Lx ~ then 

(24) D(p) = ^ y «>(p • q, q)q ® qdq. 

We shall use £ to denote either the scattering operator on the right side of eq. 1)17^ or the diffusion 
operator on the right side of eq. ()21() . The self-adjoint operator C is non-positive definite and 
represents various decoherence effects due to random fluctuations of the medium, see H27|) . The 
evolution by Eq. ()17() or H21|l preserves the mixed-state structure and hence the positivity 

(HI, (CSl). 

We shall refer to regime (a), for either H17|) or H21() . as the longitudinal case and regime (b) as 
the transverse case because the longitudinal and the transverse scales are dominant, respectively. 
Regime (c) is the borderline case. We will consider hereafter the longitudinal and transverse cases 
only. 

We note that the additional restrictions of ^ 4/3^^ and d > 3 in regime (b) are due to 
technical reasons. We believe that the results should hold for any <C Lx, d > 2 along with the 
other assumptions used in proving the above scaling limits (see [21, ^HI for details). The mean field 
result for a Gaussian potential with = and d> 3 has been established previously |22I) |Hj (see 
also '2P). 

3.2. Nonlinear kinetic equation. When U = g\'$\'^'^ , the convergence of the above scaling limits 
is not known. But because the phase of ^ is canceled in the expression of U one can reasonably 
expect that, at least before any singularity formation, U has significantly less oscillation than ^ 
and V{zLz,:x.Lx) as appearing in the NLS equation (pQ). The property of separation of scales in U 
and ^ holds in the linear case and will be assumed to hold prior to singularity formation. 

Hence we will assume below the validity of the above scaling limits for the nonlinear case and 
use the kinetic equation p7|) or (|2H) with U = g\'^\'^'^ to investigate the issue of random diffraction 
and singularity formation in the case of self- focusing (g > 0), (super)-critical {da > 2) nonlinearity 
as well as the question of rate of spread in all situations. Such a model with the kernel (|19|) was 
considered in The Fokker-Planck equation with the diffusion coefficient p2)) is an example 

of those which describe open quantum systems in contact with a heat bath of linear oscillators in 
thermal equilibrium (see, e.g., |5, 0, P, [SSI)- ^^^^ context, z would be the physical time, x 
the physical coordinates and p the momentum of the quantum particle. We will discuss the master 
equation for the open system including dissipation mechanism in the final section. 

Let us state the nonlinear kinetic equation which we will analyze subsequently: 

(25) ^' P) + P • ^^W{z, X, p) + UoWiz, X, p) = CW{z, x, p) 



where Uq is given by (jllj) and C is either the hnear scattering operator 

CW = 27tJ if(p,q)[t^(z,x,q)-T^(z,x,p)]dq 
or the hnear diffusion operator 

CW = Vp •DVpVF(z,x,p). 
The weak formulation of H25|) is given by 

(26) ^ J eVFdxdp - Jp- V^eWdxdp - J V^U • VpGVFdxdp - ^eVFdxdp = 

for smooth, rapidly decaying test functions 0. The evolution preserves the positivity of the initial 
condition which is consistent with the fact that the Wigner measures are positive and will be 
assumed subsequently (see the comment following (fT^ ). 
The nonlinear kinetic eq. (|26(1 preserves the total mass, i.e. 

-^N = 0, N = J W{z,x,p)dxdp 
but in general decreases the L^-norm 

(27) ^ J |VFp(z,x,p)(ix(ip = J WCWd:)^dp<0 

because the operator C is non-positive definite. The inequality (|27|) expresses certain irreversibility 
as a result of the weak convergence of solutions |^. One can absorb the effect of the total mass N 
into g by the obvious rescaling of W in the eq. I|26j) . Henceforth we assume that 



N = 1 

which is the case when in 

J |^„(x)|2dx=l 

for P- almost all a. 



3.3. Local existence. A natural space of initial data and solutions is the space S of the non- 
negative measures with square integrable density W 



finite Dirichlet form 



/ 



iPFpdxdy < oo, 
WCWd:!^dp < oo 



and finite, positive variances 

K = j \x- xpVFdxdp = S^- |xp G (0, oo), Vp = j \p- pl'^Wdxdp = Sp - |pp G (0, oo) 
where Sx , Sp are the second moments 

(28) S"^ = y IxpH^dxdp G (0,oo), Sp = j \p\'^Wdxdp £ {0,oo). 

In addition, we shall also assume that the initial data have a finite Hamiltonian H 
H=\Sp- ^ / p^+Ux G (-00, oo), p(x) = / Wdp 



the first term of which is the kinetic energy and the last two terms are the potential energy. A 
finite Hamiltonian and a finite variance Vp then imply 

(29) j p'^+^dx < oo, 

namely a finite potential energy. 

We recall that for the NLS eq. without a random potential and the scaling @ the local existence 
in the analogous space of functions (i.e. ^',x^', V^*, j^'l'^"''^ G L^(M'^)) holds under the additional 
assumption 
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0<cr<-^ — -, ford>3 

and a G [0, oo), for d <2, see j23j and the references therein. 

The local existence for the nonlinear kinetic equation (|26j) in the space S is assumed throughout 
the paper and will be addressed elsewhere. 

3.4. Energy law. Next, we discuss the evolution of the Hamiltonian H and the variance Vx- Let 
us first note the result of applied to the quadratic polynomials. 





(30) 


ZY-yX 


(31) 




(32) 




(33) 


U^x ■ p 


(34) 





e''^''qU{dci) = V^U 

giq x^ . qC/((iq) = X • Vxf/ 

e^'i ''2p • qf7(dq) = 2p • V^U. 

It is noteworthy that the results of the calculation are independent of 7 > and identical to those 
for 7 = (see more on this in Conclusion). 
Consider the mean dynamics for 

(35) x = j xWdxdp, p = j pWdxdp 
with the mean Hamiltonian defined as 

(36) ^ = ^|Pl'- 

Using the above and integrating by parts we obtain the following 



d 

—p = I C*pWdxdp 
oz 

d 



'xdp + j . 



^Yp = j V^U -pWdxdp^ I C*\p\^Wdpdx 
p'^+^dx = [ VJJ ■ pWdpdx 



d g 
dza + l 
and hence 

(37) = p - [ C*pWdxdp 

oz J 

(38) = J C*\p\'^Wdpdx 

which will take a more explicit form once we calculate C*p, >C*|pp with C* of each case. 



As is the case without a random potential and the scahng hmit Q) 23j, we assume subsequently 
that the laws ()37() . (|38() for the Hamiltonian holds for the local solutions of eq. ()26() in the space S. 

In the following sections, we first derive the variance identity for eq. H26() in the longitudinal 
case (regime (a)), then the transverse case (regime (b)) and finally the case with friction. As is the 
case without a random potential and the scaling limit (jlj) [^, we assume throughout the paper 
that the variance identity derived below holds true for the maximally extended local solutions of 
eq. (|26|) in the space S. Let [0, z*) be the interval on which the maximally extended local solution 
is defined. When = oo then the local solutions become global solutions; when z^, < oo then 
the solutions are said to develop finite time singularity. Our main goal then is to analyze what 
these formally derived relations tell us about singularity formation and rate of spread of low Fresnel 
number nonlinear parabolic waves in random media. 

Before concluding this section let us state some elementary inequalities which will be useful later. 

An application of the Cauchy-Schwartz inequality and the marginal positivity 1)141) . (|15() to the 
first moments x, p then leads to 

(39) |x|2 < Sx, |p|2 < Sp. 

Furthermore, by using the mixed-state structure H13|) in estimating the cross moment 



(40) Sxp = y X • pWdxdp. 
one deduces that 

(41) Sip < SxSp < cx). 
Likewise the covariance 

Vxp = y (x - x) • (p - p)Wdxdp = Sxp-x-p 

can be bounded by Vx and Vp as 

(42) < VxVp. 

4. The longitudinal case 

4.1. Variance identity. The variance identity has been long used to derive the wave collapse 
condition for the NLS without the random potential . Below we reformulate it and extend it to 
the phase-space models (j^ . 

We have the following simple calculations for C* = C: 

Cx = 

Cp = 27ry"$(0,p-q)(p-q)dq = (hv (|TH|)) 



C\x 
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£x-p = x-27ry$(0,p-q)(p-q)dq = fbv ^) 

£|p|2 = -27ry"$(0,p-q)[|p|2-|q|2]dq 

= -27r j $(0,q)(2p-q)-qdq 

= 27r y $(0,q)|qpdq = R. 



We have used the hnear Bohzmann operator L in the above calculation; the same result holds for 
the linear diffusion operator £ for which case 



i? = 2 X trace(D) 

where the diffusion matrix D is given by H22|) . We shall use the above identities to perform 
integrating by parts in the derivation of the variance identity. 

The evolution of the mean position x and momentum p is then given by 

d 
oz 

= 0, 

as a consequence the mean Hamiltonian H is invariant. Moreover, the evolution of (x, p) for the 
autonomous Hamiltonian system of the harmonic oscillator. 
The evolution of the variance Vx is given by 

d 



Differentiating Sxp we obtain 



—Sxp = Sp- -^-^ j p + dxdp 
Hence the second derivative of Vx becomes 

(43) = 4,^,-5) + ?(^/,-.x. 

An alternative expression for the variance identity is 

(44) ^Vx = 2da{H -H) + {2-da)Vp. 

Both forms (|43j) and (|44|) will be used to obtain dispersion estimates below. 

4.2. Dispersion rate. Although the medium is lossless, reflected in the fact that the total mass 
= 1 is conserved, but the Hamiltonian is not conserved by the evolution since the scattering 
with the random potential is not elastic. Indeed, its rate of change is 

d 

(45) —H = R, H{z) = F(0) + Rz 

due to the diffusion-like spread in the momentum p. 
In the critical case dcr = 2, we have the exact result 

^Vx = m- 2|p|2 = A{H{Q) -H)+ 4Rz 

before any singularity formation and hence the following. 
Proposition 1. If da = 2 or g = 0, then 

(46) Vx{z) = Vx{0) + 2Vxp{0)z + 2{H{0) - H)z^ + ^z\ z G [0, z,). 
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The analogous result {Vx ~ z"^) for the linear Schrodinger equation {d = 1, g = 0) with a random 
potential has been proved previously [7j. 
We have from that 

and hence 

—rVx < iH -4H = 4:(H(0) + Rz) - AH, for g(2 - da) < 
oz^ 

— K > 4:H-AH = 4.(H(0) + Rz) - AH, for g(2 - da) > 0. 
On the other hand, from (|44() we obtain for any g 

— K < 2da(H-H), for 2 - da <0 
oz^ 

— > 2da(H - H), for 2 - da > 0. 

oz^ 

Integrating the above inequalities twice, we obtain the following. 
Proposition 2. The estimates hold 

Vx{z) < T4(0) + 214p(0)z + 2{H{{)) - H)z'^ + '^Rz'^ , for g{2 - da) < 



Vx{z) > Vx{0) + 2Vxp{0)z + 2{H{0)-H)z^ + ^Rz\ forg{2-da)>0 



and 



Vxiz) < Vx{0) + 2Vxp{0)z + da{H{0)- H)z'^ + '^Rz^, for 2 < da 

da 



Vx{z) > VxiO) + 2VxpiO)z + da{H{0)-H)z^ + —Rz^, for 2 > da 
for all z £ [0,z^). 
Therefore 

Corollary 1. Assume g <0 (hence H >0). Then 

Vx{z) < K(0) + 2Kp(0)z + {da V 2)[H{<d) - H]z^ + ^^i^z^ 

and 

(48) Vxiz) > T4(0) + 2VxpiO)z + (da A 2)[H{0) - H]z^ + ^^i^z^ 
for z £ [0,z^) 

4.3. Singularity formation. Finite-time singularity for the critical or supercritical NLS {da > 2) 
without a random potential is a well known effect In this case the singularity is the blow-up 
type Vp, 00. Here we take as a model equation to gain some insight into singularity 

formation in the presence of a random potential. 

First we consider the self-focusing case 5 > 0. For da > 2 one can bound as 

(49) Vx{z) < VxiO) + 2Vxp{0)z + da{H{0) - H)z^ + ^z^ = F{z) 
and looks for the situation when F{z) vanishes. 
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The sufficient conditions for F{z) to vanish at a finite positive z are that F{z) takes a non- 
positive value F{zo) < at its local minimum point zq > 0. The local minimum point zq is given 
by 



H - H{0) + \ iH{0) - HY - 2RV,p{0)/{da) 
(50) ZQ = 



R 

Therefore we are led to the singularity conditions for g > 0. 

Proposition 3. For da > 2,g > 0, the local solution of the nonlinear kinetic eg. ii^6|) with 
the longitudinal randomness develops singularity at a finite time < zq given by H5(J\) under the 
condition F{zo) < and either one of the following conditions 

(51) V,p{0) < 

(52) V^piO) > 0, H{0) < H - 



2RV,p{0) 



da 

Remark 1. Clearly, the condition F{zq) < requires H{Q) to be sufficiently below H by allowing 
the potential energy 

^ ( p^+^dxdp 



a+1 

at z = to be sufficiently negative. We will give an explicit blow-up condition for the supercritical 
nonlinearity in the next section. 

In the linear or self-defocusing case g < the right side of or is always positive in view 
of the inequality 

2|Kp| <V:, + Vp, 

cf 



Therefore, the above result suggests that the scattering term does not prevent large scale singu- 
larity, even though it appears (from the upper bound zq) that it takes longer time than the time 
to singularity in the homogeneous NLS |23j . 

With the assumption that the variance identity and the energy laws (|37|) . (|38|) hold for the 
maximally extended local solutions of eq. (|26j) then it is clear that the singularity is the blow-up 
type 



or equivalently 



lim / p^^^d:>c = oo 



lim Sp = oo 

Z—*Zf 



by the finitude of the Hamiltonian. 

Next we will follow the argument of |14j to show more explicitly the blow-up mechanism in the 
case with the supercritical, self-focusing nonlinearity and give a sharper bound on z* under certain 
circumstances. 

Proposition 4. Suppose da > 2, g > 0. Then under the conditions 115 
(53) H{zo) = H{0) + Rzo<0 

with Zq given by i5(jl\) . as well as Vp blow up at a finite time where 

2K(0) 



< Zq A 



V^p{0){2 - da)- 
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Let us give the argument below. Since blow-up is a local phenomenon, is a poor indicator of 
its occurrence. To this end a more useful object to consider is Vp. 
From (gll) it follows that 

(54) -^V^ < {2-da)Vp<0, z < z,. 

Hence V^{z) is a negative, decreasing function for z < z*. Also by the Cauchy-Schwartz inequality, 
we get 

(55) 0<V^p< V^Vp < V,iO)Vp, 0<z<z, 
and hence 

Let A{z) = —V^{z) > 0,2 < z*. We have from (|44j) and (|55|) the differential inequality 



which yields the estimate 



^(0) 1 



1 - zCA{0) ' CA{0) 

and thus the blow-up of —V^{z) at a finite time. This along with H56|) then H53j) implies the 
divergence of at a finite time 

lim Vp(z) = oo 

with 

1 2T4(0) 



(58) < 



C7^(0) V,pi0)i2-da)' 



For high power dcr ^ 2, ((S5|) is a better upper bound for z^ than zq given by (I^U)) . 

The preceding argument demonstrates clearly the blow-up mechanism, namely the quadratic 
growth property ((S7|l as well as gives a sharper bound on the blow-up time for large da. 



5. The transverse case 

The scattering operator C in this case corresponds to elastic scattering, instead of the inelastic 
scattering of the longitudinal case. This affects the evaluation of £* = £ on quadratic polynomials: 

£x = 

|2 |^|2 



£p = 2^y5(^Li^)y"$(^i;,p-q)du;(p-q)dq = 

£|x|2 = 
£x • p = 



£|p|2 = -27r J 6{ 'P'' ^ '^^l' ) J ^w, p - q)dw [|p|^ - |q|^J dq = 0. 
The same results hold for the diffusion operator £ with the diffusion matrix H23|) . 
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5.1. Rate of dispersion. The variance identity still holds Hence the second derivative of Vx 
becomes 

(60) = 2da{H - H) + {2 - da)Vp. 

The main difference from the longitudinal case is that the Hamiltonian is invariant 

oz 

By the same argument as in the longitudinal case we have the following analogous estimates. 
Proposition 5. If da = 2 or g = 0, then 

(61) Vx{z) = 14.(0) + 2Vxp{0)z + 2{H{0) - H)z\ z E [0, z.). 
Proposition 6. The following estimates hold for z E [0, z^): 

Vx{z) > K.(0) + 2T4p(0)z + 2{H - H)z^, for g{2 - da) > 
Vx{z) < Vx{0) + 2VxpiO)z + 2{H - H)z\ for g{2 - da) < 

and 

Vx{z) > VxiO) + 2VxpiO)z + da{H - H)z^, for 2 > da 
Vxiz) < VxiO) + 2VxpiO)z + da{H - H)z'^, for 2 < da. 

Therefore, 

Corollary 2. Assume g < (hence H >0). Then for z E [0, z^,) 

(62f^.(0) + 2Vxp{0)z + {da A2)[H - H]z'^ < V^iz) < K(0) + 2Kp(0)z + {da V 2)[H - H]z^ . 
From these estimates we see that a ballistic kind of motion takes place in the transverse case. 

5.2. Singularity. By finding the zeros of the upper bounds in Proposition |H1 we can derive the 
conditions for singularity formation. As in the longitudinal case, for 5 < 0, 1)61(1 and the left side 
of (|62|) are always positive. 

Proposition 7. For g > 0,da > 2, the local solution of the nonlinear kinetic eq. i2b}) with the 
transverse randomness develops singularity in a finite time z^, < 00 under either of the following 
conditions 

(63) H <H 
or 

(64) H>H, Vxp{0)<0, H <H+\Vxp{0)\^/{daVx{0)). 
Remark 2. It is easy to see that the singularity time z^ is bounded from above by 



^''> ^° = MH-ii) 

under condition i6cl\) and by 



-Vxp{0) - ^|Kp(0)|2 - daVx{0){H - H) 
" da{H - H) 

under condition \64^ - 

14 



The preceding singularity conditions are identical to those for the homogeneous (super)-critical 
NLS equation |2S1- This is already suggested by the previous result where is shown to 
be the instability condition for the diffusion approximation of the kinetic equation (|26() with the 
kernel (|T^. 

With more stringent conditions one can demonstrate a stronger sense of blow-up. Let us state 
the result. 

Proposition 8. Suppose da > 2, g > 0. Then under the conditions 

H <H, 

as well as Vp blow up at a finite time z^, < zq. 
Let us sketch the argument below. From H6U|) it follows that 

(66) ^14. < {2-da)Vp<0, z<z,. 

Since V^{z) is a decreasing function and becomes negative after a finite time z when Vx{z) achieves 
its maximum V^. 

Again for A{z) = — V^'(z) > 0, z > 2 we have from (|66() and (|55() the differential inequality 

9.^„.2 ^ da — 2 
—A > CA^ C = —— > 0, z>z. 
oz 4:Vx{z) 

This is a plausible mechanism for the finite-time blow-up of the kinetic energy and the potential 
energy. 

6. Fokker-Planck equation with dissipation 

In this section we will apply the same analysis to the Fokker-Planck equation (|26|) with, instead, 
the Caldeira-Leggett operator |1] 

(67) CW = DApW + XVpi-pW], D,X>0 

which is, perhaps, the simplest example of the phase-space models for open, dissipative systems 
[S], pp. As remarked before in this context z is the physical time, x the physical coordinates and p 
the momentum in R^. If the open quantum system is the electron motion in a semi-conductor then 
the nonlinearity arises from the electrostatic Hartree potential self-induced by the density p with 
self-defocusing nonlinearity. Since operator 1)6 7|) is not self-adjoint (w.r.t. the Lebesgue measure) 
it is outside the scope of the NLS eq. with a real- valued random potential (^. The method used 
here can be applied to such a problem. However, in order to compare how different decoherence 
mechanisms affect wave spreading and collapse we will still consider the same power-law nonlinearity 
as in the previous cases. 

The dissipative effect is of frictional nature and the Fokker-Planck equation with H67() still pre- 
serves the total mass 

oz 

As before we set = 1. 

We re-calculate the action of £*, the adjoint operator, on the quadratic polynomials: 



£*x = 

£*p = -Ap 

£*|x|2 = 

C*x • p = — Ax • p 

C*\p\^ = 2dD-2X\p 
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6.1. Rate of Spread. Using these results and integrating by parts we get 

d__ 

d__ 

which again constitute an exactly solvable system: 



Ap 



(68) 
(69) 



P(^) 



-\z 



P(0) 

x(0) + PM(i_e-^-). 
A 



The Hamiltonian H, defined by H36|). of the mean motion is then given by 

H{z) = H{{))e-'^^'. 

Also, 



oz 
dz 



2K 



xp 



da , 



da{H -H) + {1- —)Vp - AKp 



= 2{H-H) + 
From these we obtain the upper bound 

(70) V.,p{z) < V^p{Q) + dae-^' 

(71) V.,p{z) > V^p{0) + dae-^' 

and the bounds 
(72) 

(73) 



(2 - da)g /" +i 



a + 1 



V^p{z) < V^p{0) + 2e 



e'^{H{s) - H{s))ds, for da > 2 
e'^{H{s) - H{s))ds, for da < 2 

e'^{H{s) - H{s))ds, for ^(2 - da) < 
e'^{H{s) - H{s))ds, for 5(2 - dfj) > 0. 



On the other hand, the rate of change of the Hamiltonian is given by 
d 



dz 



[H{z) - H{z)\ 



dD - AK 



dD - 2\{H - H) 



2Xg_ 

a+1 

For self-focusing nonlinearity g > the upper bound on H follows 



p^+^dx. 



(74) 



H{z) - H{z) < H{0) - H{0) + (1 - e 



-2A2 



dD 



2A 



whereas for the self-defocusing nonlinearity g < we have both the lower as well the upper bounds 



H{0) - H{0) + (1 - e-^^')— < H{z) - H{z) < H{0) - H{0) + dDz. 

2a 



dD 



Therefore for g >0 we have the upper bounds 
V^piz) < Kp(0) + da 



g(0)-g(0) ^^ 
A 
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-Xz 



dD 



whereas for g <0, we have the esthnates 
V^p{z) > K.p(0) + (da A 2) 

V^p{z) < Kp(0) + (da V 2) 

Therefore, 

Proposition 9. For g > 0, z [0, z*) 



A 

H{0)-H{0) dD 
A 



(1 



' 2A2 

, cfaD 
) + ^z. 



(75) 



da 



dD^ 



Vxiz) < V,{0) + 2V,p{0)z + —{2H{0)-2H{0) + —)z 



2da 



A2 



A 

H{0) - H{0) + ^ 



(l_e-A-) + 



A 

2A3 ^ 



whereas for g < 0, z £ [0, z^ 

;(0)+2 
2((ia A 2) 



VAz) > V,{0)+2V,p{0)z + ^^{2H{0)-2H{0) + ^)z 



A2 



^(0) - H{0) + ^ 



1 



-A.) I (daA2)dZJ ^^ 



-2A2^ 



2A3 



^x(^) < F.(0)+2y,„(0)z + 2(daV2 



ri/(o) 


dD' 




A 







-Az 
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+ 



d^aD 



A 



6.2. Singularity. The singularity conditions for the self-focusing case are determined from that 
the right side of (|75jl has a positive root. 

Proposition 10. For da > 2, g > 0, the local solution of the Fokker-Planck eq. h2b\) with \61\) 
develops finite-time singularity if 



(76) 



l^xp(0) + y(/f(0)-^(0) + ^)<0. 



The condition (jTSJ is better suited for large A and can be much improved for small A. As A ^ 
this case should approach to that of the longitudinal case with D given by (|^^ being a scalar. 
Let us sketch an explicit blow-up mechanism in the super-critical case under different conditions. 

Proposition 11. Suppose da > 2, g > 0. Then under the conditions 



(77) 
and 
(78) 



V,p{0) < 0, HiO) - 5(0) + ^ < 



K(0) < 



da 



2A 



as well as Vp blow up at a finite time z^, where 



z^ < — In ■ 



da 



A do--2-2Ay^(0)' 
We present the argument here. Rewriting some of the preceding calculations, we have 

d 



and 
(79) 



= da{H -H) + {1- ^)Vp - XV^p. 
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da , 
Y 



Condition ((TTj) and ^ imply that 

Using this in ((7^1) we obtain 

d_ 

dz 

which together with ((77|) impUes Vxp{z) < 0,z < z* and Vx[z) is decreasing. Using the inequahty 



H{z) - H{z) < -e^''^^. 



Vr,> 



" " K(o) 

in H79|) we obtain again the quadratic growth estimate 



for ^ = — T4p > 0, or equivalently 



for 

From this it follows that 



9 „ da — 2 \,^o 



B = e^'A. 



1 da-2 p-\z\ 

the right-hand side of which blows up at a finite time if and only if (|78j) holds. 

The finite-time blow-up of B would also lead to the finite-time blow-up of A and hence of Vp. 

7. Conclusion 

We have derived several kinetic equations meant to model the large scale, low Presnel number 
behavior of the NLS equation with a rapidly fluctuating random potential based on the rigorous 
theory [S], jlUj for the linear Schrodinger waves in the same situation. The low Fresnel number 
waves interacting with a rapidly oscillating potential give rise to, in the scaling @, a self-averaging 
limit of a deterministic kinetic equation with a scattering operator. 

We analyze these kinetic equations in order to shed light on two problems: the rate of dispersion 
and the singularity formation. Our main assumptions are the existence of local solutions of the 
kinetic models in the space 5, that the local solutions satisfy the Hamiltonian laws (jSZl)) jHEl) and 
the variance identity. 

The scattering operator in the kinetic equations may be a Boltzmann-type operator or a Fokker- 
Planck operator, depending on whether the random potential is fully resolved by the wave or not. 
This, however, does not affect our results for either problem. 

What is more important for our investigation is the structure of the scattering kernel in the 
Boltzmann operator or the diffusion matrix of the Fokker-Planck operator. There are two types of 
structures: the longitudinal regime when the longitudinal scale of the random potential is dominant 
and the transverse regime when the transverse scale of the random potential is dominant. The 
third kind of kinetic equation that we have analyzed is based on the Fokker-Planck equation with 
the longitudinal type of diffusion matrix plus a frictional term. This model is motivated by the 
Caldeira-Leggett operator for open quantum dynamics. 

For the problem of dispersion, we have shown that the kinetic equations of the longitudinal type 
produce the cubic (z^) law, that the kinetic equations of the transverse type produce the quadratic 
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(z^) law and that the Fokker-Planck equation with friction produces the hnear (z) law for the 
variance Vx{z) for z S [0,2:*). 

For the problem of singularity, we have shown by analyzing the variance identity that the singu- 
larity and blow-up conditions in the transverse case remain the same as those for the homogeneous 
NLS equation with critical or supercritical self-focusing nonlinearity, that the singularity and blow- 
up conditions have changed in the longitudinal case and in the frictional case. 

Finally let us make a slight extension of our results by considering the nonlinear kinetic equation 
of the more general form 

(80) ^' P) + P • VxVF(2, X, p) + U^W{z, X, p) = CW{z, X, p) 

with 7 > small, but positive and the mixed-state Wigner distributions H13|l as initial data. Eq. 
(|26)) is formally the geometrical optics limit of eq. (|80)) . Eq. (|8fl|) retains the diffraction effect of 
the nonlinear potential term that is absent in eq. (|26)) with 7 = 0. Even though eq. (|80j) does 
not preserve positivity, it can be shown to preserve the mixed state structure (|13)) and therefore 
the marginal positivity (|14j) . (|15j) as well as the inequalities (|41|) . (|42j) in both the longitudinal 
and transverse cases. Since the calculations (|30j) - (|34|) are independent of 7 > 0, all the results 
in Section 4 and 5 hold for eq. (|80)) with 7 > as well. However, because the Caldeira-Leggett 
operator 1)6 7|) is in the Lindblad form only in the classical limit 6^, therefore in order to maintain 
the solution in the the mixed-state form (jl^j) we need to set 7 = in this particular case. 
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